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ABSTRACT 

We present the first detailed study of the stellar populations of star- forming galaxies at z ~ 1.5, 
which are selected by their [O n] emission line, detected in narrow-band surveys. We identified ~1,300 
[O n] emitters at z = 1.47 and z = 1.62 in the Subaru Deep Field with rest-frame EWs above 13A. 
Optical and near- infrared spectroscopic observations for «10% of our samples show that our separation 
of [O n] from [O in] emission-line galaxies in two-color space is 99% successful. We analyze the multi- 
wavelength properties of a subset of ^1,200 galaxies with the best photometry. They have average 
rest-frame EW of 45A, stellar mass of 3 x 10 9 M Q , and stellar age of 100 Myr. In addition, our SED 
fitting and broad-band colors indicate that [O n] emitters span the full range of galaxy populations at 
z ~ 1.5. We also find that 80% of [O n] emitters are also photometrically classified as "BX/BM" (UV) 
galaxies and/or the star-forming "BzK" (near-IR) galaxies. Our [O Ii] emission line survey produces 
a far more complete, and somewhat deeper sample of z ~ 1.5 galaxies than either the BX/BM or 
sBzK selection alone. We constructed average SEDs and find that higher [O n] EW galaxies have 
somewhat bluer continua. SED model-fitting shows that they have on average half the stellar mass 
of galaxies with lower [O n] EW. The observed [O n] luminosity is well-correlated with the far-UV 
continuum with a logarithmic slope of 0.89 ± 0.22. The scatter of the [O n] luminosity against the 
far-UV continuum suggests that [O n] can be used as a SFR indicator with a reliability of 0.23 dex. 

Subject headings: galaxies: photometry — galaxies: distances and redshifts — galaxies: evolution — 
galaxies: high-redshift — infrared: galaxies — ultraviolet: galaxies 



1. INTRODUCTION 

Since the star formation rate (SFR) density of galax- 
ies is an order of magnitude higher at z w 1-5 than 
in the local universe (e.g., Hopkins 2004), emission- line 
galaxies are prominent in the young universe, and thus 
are useful probes of the evolution of galaxies. Current 
techniques to identify such emission-line galaxies include 
grism surveys (Shim et al. 2009; Atek et al. 2010; Coil 
et al. 2011; Straughn et al. 2011), slit spectroscopy (e.g., 
DEEP2; Davis et al. 2003), and narrow-band imaging 
(e.g., Malkan et al. 1996; Tcplitz et al. 1998; Fujita et 
al. 2003; Hippelcin et al. 2003; Ly et al. 2007, here- 
after L07). Narrow-band surveys have identified large 
samples of star-forming galaxies by detecting their red- 
shifted nebular emission lines, at low and intermediate 
redshifts, to as high as z ~ 7. These surveys have the 
ability to: (1) determine rcdshift to 1 percent accuracy, 
(2) derive emission-line fluxes, which can be used to de- 
termine star-formation rates (SFRs) and trace the evo- 
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lution of the cosmic SFR density. And they accomplish 
this with roughly an order-of-magnitude higher efficiency 
than spectroscopic surveys. 

However, it is still unclear what galaxy population(s) 
these techniques probe or are unable to probe. For ex- 
ample, it might be that these surveys are biased towards 
the youngest and possibly least reddened galaxies. Such 
claims have been made for Lya emitters when compared 
to Lyman break galaxies at similar redshifts (e.g., Ga- 
wiser et al. 2007). Thus, it is crucial to understand the 
selection bias of current and future emission-line surveys 
using these techniques. Our survey utilizes ~1,300 [O n] 
emitters at z = 1.47 and z = 1.62 in the Subaru Deep 
Field (SDF; Kashikawa et al. 2004). 

The outline of this paper is as follows. In Section 2, 
we discuss the SDF and the optical imaging that we have 
acquired and reduced. Section 3 discusses the selection 
of narrow-band (NB) excess emitters, the spectroscopic 
observations, the technique to separate [O n] emitters 
from Ha and [O III], and the photometric redshifts. Sec- 
tion 4 presents the results on the [O n] EW distributions, 
a comparison between the rest-frame UV luminosity and 
the [O n] luminosity, their stellar population from model- 
ing the spectral energy distribution (SED), and compares 
the [O n] population against popular color selection tech- 
niques at these redshifts. We then summarize our con- 
clusions in Section 5. Throughout this paper, we adopt 
a flat cosmology with = 0.7, CIm = 0.3, and H = 70 
km s~ 4 Mpc" 1 , and magnitudes are reported on the AB 
system (Oke 1974). 

2. THE SUBARU DEEP FIELD OBSERVATIONS 
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Several extragalactic fields already have deep imaging 
in NB filters with Subaru/Suprime-Cam (Miyazaki et al. 
2002). Of all of these, the SDF has the most sensitive 
optical NB imaging in six NB filters in the sky, and is 
further complemented with deep multi-band imaging be- 
tween 1500A and 4.5^m. The most prominent emission 
lines entering these NB filters arc Ha, [O in] A5007, and 
[O n] AA3726,3729, at well-defined redshift windows be- 
tween z ~ 0.07 and z m 1.62. 

2.1. NB921 and NB973 Imaging Observations 

For this paper, we focus on excess emitters found us- 
ing two filters, NB921 and NB973. Their central wave- 
lengths and FWHM are 9196A and 132A for NB921 and 
9755A and 200A for NB973. The acquisition of the SDF 
NB921 data has already been discussed in Kashikawa et 
al. (2004), Kashikawa et al. (2006), and L07. The NB921 
data are identical to the publicly released SDF v.l data. 10 
The NB973 data were acquired on 2005 March 16-17 and 
2007 May 9 for a total integration time of 21.5 h (lye et 
al. 2006; Ota et al. 2008). The seeing varied between 0'.'5 
and l'.'O and observing conditions were photometric in 
2005 but not in 2007. 

Both the NB921 and NB973 data were reduced using 
SDFRED (Yagi et al. 2002; Ouchi et al. 2004), a software 
package designed solely for Suprime-Cam data. The spa- 
tial resolutions for the final NB921 and NB973 mosaics 
are 0'.'96 and 0'.'91, respectively. 

By repeatedly measuring the flux within 2" diameter 
apertures at random positions in the mosaics, we deter- 
mined that the 3cr limits for the NB921 and NB973 data 
are 26.60 and 25.64 mag, respectively. A summary of the 
sensitivity in other wave-bands, that are later discussed, 
can be found in Table 1 of Ly et al. (2011). The pub- 
licly available z'-band mosaicked image, which is used 
throughout this paper, has a 3c sensitivity of 26.83 mag, 
while the z'-band data (later used) is 0.93 mag deeper. 

3. SELECTION OF NARROW-BAND EXCESS EMITTERS 

The NB921- and NB973-detected catalogs were ob- 
tained by running SExtractor (Bertin & Arnouts 1996, 
vers. 2.5.0). The unmasked regions cover 870.4 arcmin 2 
for NB921 and 788.7 arcmin 2 for NB973. The latter 
is smaller due to higher systematic noise in one of the 
ten CCDs, which was masked to avoid significant spu- 
rious detections. A total of 119,541 and 84,786 sources 
were detected in the unmasked regions of the NB921 and 
NB973 images, respectively. 

To select NB excess emitters, the standard technique 
is to compare the flux in the NB against the flux in an 
adjacent broad-band filter. We use the z'-band, since 
it samples the local continuum around these NB filters. 
Fluxes from the z'-band were extracted by running SEx- 
tractor in "dual-image" mode, where the respective NB 
image was used as the "detection" image. This works 
well because all three mosaicked images have very simi- 
lar seeing, so the z'-NB colors are determined within the 
same physical scale of the galaxies. For the extraction 
of fluxes and selection of sources, we use a 2" diameter 
circular aperture, which encloses 17 kpc at z w 1.5. 

10 http: //soaps .nao . ac . jp/SDF/vl/index . html. 
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Fig. 1.— Total system throughput for the z' (solid), NB921 
(dashed), and NB973 (dotted) filters. 

One minor issue affecting the NB921 and NB973 filters 
is that they are located redward of the z'-band sensitiv- 
ity peak, ~8700A (see Figure 1). For redder galaxies, 
this makes the z'-NB colors systematically higher, and 
could falsely identify lineless red galaxies as NB-excess 
emitters. 

To resolve this issue, we empirically estimate the con- 
tinuum flux at the wavelengths of the NB filters by fitting 
the i'—z' and z'-NB colors for non-NB excess sources: 



(z'-NB921)cr 
(z' - NB973) C : 



-0.003 
-0.014 



0.0770'' 
0.271(i' 



z'), and 
z>). 



(1) 
(2) 



These relations are consistent with simple linear extrap- 



olations of the 



color to 9200A or 9755A. The cor- 



rected z'-NB color is then 

A(z' - NB) = (z' - NB) - (z' - NB) C . 



(3) 



To illustrate the reliability of these corrections, we plot 
the z'-NB973 colors with and without the color-term cor- 
rection for NB973 emitters in Figure 2. This correction 
was not performed previously for the NB921 excess emit- 
ter sample identified by L07. It has a minor effect for 
the NB921 sample, but is more important for the NB973 
sample since the NB973 filter is significantly redder by an 
additional ~560A. For example, we find that the rms in 
the z'-NB color (at the bright end) is decreased by 10% 
for NB921 and 50% for NB973 when this correction is 
applied (see Figure 2). Since the minimum threshold for 
selection (discussed below) is driven by the observed scat- 
ter, such a reduction permits us to select toward lower 
minimum z'-NB excesses. 

In Figure 3, we illustrate the corrected z'-NB921 
and z'-NB973 colors for NB921- and NB973-detected 
sources, respectively. For selecting objects with genuine 
NB921 and NB973 excesses, we require minimum z'- 
NB magnitude differences of 0.20 mag (observed EW of 
32 A) and 0.25 mag (observed EW of 70 A), respectively. 
These minima are 5.6 and 6.7 times the rms of the z'-NB 
color at the luminous end. The observed EW is deter- 
mined from a combination of the broad- and narrow-band 
fluxes, and is given by: 



EW obs = ANB 



A,NB 



h 



L/a,,-/a,nb(ANB/Az) 



(4) 
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Fig. 2. — z'— NB973 colors for NB973 excess emitters as measured (left) and after applying the i'—z' color-term correction given in 
Equation (2) (right). Cyan, black, and red squares denote galaxies with blue, intermediate, and red i'—z' colors. The color-term correction 
ensures that red galaxies with weak or absent [O II] emission are not falsely classified as NB973 excess emitters. 




Fig. 3. — A(z'-NB) for NB921 (left) and NB973 (right) samples (grey points) in the SDF. The 3a and 5a sensitivity limits for the NB data 
are shown as vertical dashed lines. Minimum excesses of 0.2 mag (left) and 0.25 mag (right) arc illustrated with the horizontal dashed lines. 
Also overlaid are curves for TV X o"^/_ NB with N = 2.5 (dashed line) and 3.0, 4.0, and 5.0 (solid lines). Sources that fall above the diagonal 

dotted lines are undetected in the ,z'-band at 3<r. In the absence of noise, the strongest possible magnitude excess is ^2.5 magnitudes, 
which corresponds to a pure emission line centered in the NB filter, and no continuum. Spcctroscopically confirmed sources are overlaid 
with symbols representing different instruments: DEIMOS (circles), Hectospec (squares), FOCAS (triangles), and FMOS (diamonds). Red 
crosses represent galaxies that are dual excess emitters in both the NB704 and NB921 bands (see Section 3.1.2). Ha, [O III], [O II], and 
Lya are distinguished by red, green, blue, and purple colors, respectively. Also overlaid in grey symbols are the cases where an emission 
line was not detected in the spectra. 
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where ANB and Az (955A) arc the FWHMs for the 
narrow-band and z', respectively, and f\x is the flux 
density in erg s _1 cm -2 A -1 in band "AT" (L07). We 
also require that the z'-NB color exceeds 



3cr 



-NB 



(NB) = -2.5 log 1 



/ia3<tNB + / 



V^OZ 



NB 



(5) 

where f Vt saX is the 3cr flux density limit for band "A"" 
and /j/,nb is the flux density for a given NB magnitude. 
These selection criteria yield 2512 NB921 excess emitters 
and 1666 NB973 excess emitters. The 3cr emission-line 
flux limits are 5.8 x 10 -18 and 1.7 x 10~ 17 erg s _1 cm~ 2 
for NB921 and NB973, respectively. Smaller samples 
were also constructed by requiring stronger detections 
of 4ct / and 5a i These samples contain 1912 

z — NB z — NB 1 

and 1616 for NB921 and 839 and 565 for NB973. 

While the color-term corrections (see Equations (1)- 
(2)) should apply to galaxies with SEDs that behave sim- 
ply, at z ~ 1.3-1.4 the 4000A break falls in the z'-band, 
but is blueward of the NB filter. This poses a problem, 
as galaxies with older stellar populations (>1 Gyr) will 
have a strong break, and the continuum flux at ~9200A 
or 9700A can easily be underestimated from a linear ex- 
trapolation of the i'—z' color. This population was con- 
firmed when the photometric rcdshifts (hereafter photo- 
z's or Zphot; see Section 3.3 for further discussion) were 
derived. At these z p hot's, the [O n] emission line does not 
fall into the NB973 passband. Since these galaxies are 
very red, simply imposing a cut of i'—z! < 0.8 (0.4) for 
NB921 (NB973), most of these 4000A-break galaxies can 
be removed without affecting the true line emitters. By 
removing these red false NB excess emitters, the NB921 
sample is reduced to 2361 (3a), 1787 (4a), and 1505 (5a). 
Likewise, the NB973 sample is reduced by about 25%- 
30%, to 1243 (3cr), 571 (4a), and 375 (5a). 

3.1. Independent Confirmation of NB Excesses 
3.1.1. Follow-up Spectroscopy of NB Excess Candidates 

To confirm that the NB technique efficiently identifies 
emission-line galaxies, we have conducted spectroscopy 
with several instruments (both optical and near-infrared) 
over the past several years. The observations were ob- 
tained with Keek's Deep Imaging Multi-Object Spec- 
trograph (DEIMOS; Fabcr et al. 2003) in 2004, 2008, 
and 2009, Subaru's Faint Object Camera and Spectro- 
graph (FOCAS; Kashikawa et al. 2002) in 2004 and 
2007, MMT's Hectospec (Fabricant et al. 2005) in 2008, 
and Subaru's Fiber Multi-object Spectrograph (FMOS; 
Kimura et al. 2010) in 2012. 

The DEIMOS and FOCAS follow-up spectroscopy has 
been discussed in Kashikawa et al. (2006) and L07, and 
DEIMOS data that were acquired more recently are dis- 
cussed in Kashikawa et al. (2011). The Hectospec and 
FOCAS observations typically had 1-2 hrs of on-source 
integration, while the DEIMOS observations had 2-4 hrs. 
In addition, we observed for 2.5-3.0 hrs with FMOS. 

For the Hectospec observations, we utilized the 270 
mm -1 grating blazed at 5200A to yield a spectral cover- 
age of 3650A-9200A with a spectral resolution of «5A. A 
total of four different fiber configurations were used. The 



on-source integration varied between 80 and 140 minutes. 
The MMT spectra were reduced following standard pro- 
cedures with the hsred 11 reduction pipeline. 

FMOS observations, which only targeted galaxies sus- 
pected to be [O il] emitters (see Section 3.2), were ob- 
tained on 2012 8-9 Apr utilizing "cross-beam switching" 
mode. Each exposure was 15 m in length and obtained 
with "ABAB" dithering. The data were reduced using 
fibre, a reduction package that consists of iraf scripts 
and CFITSIO routines, designed solely for optimal FMOS 
data reduction (Iwamuro et al. 2011). 

A summary of the spectroscopic observations obtained 
for NB921 and NB973 emitters is provided in Table 1. 
In total, 522 NB921 emitters and 204 NB973 emitters 
were targeted. Our NB973 spectroscopic observations 
are sparser, since that sample became available much 
later. The majority (70% and 90% for the NB921 and 
NB973 samples, respectively) of our spectra are from 
Hectospec and FMOS, with most of the [O n] confirma- 
tions from FMOS. Among these spectroscopically tar- 
geted NB921 excess emitters, «57% (295) have reliable 
rcdshifts (based at least one other emission line). Simi- 
larly, 51% (104) of the observed NB973 excess emitters 
yielded clear spectroscopic identifications and redshifts. 
These success rates were closer to 70% using only the op- 
tical spectroscopy, which is generally more sensitive. In 
addition, we exclude 31 (6%) NB921 and 8 (3%) NB973 
sources that show single weak emission lines in their opti- 
cal spectra which, if they are real, are consistent with the 
redshifts predicted by an emission line in the NB filter. 

The numbers of spectroscopically confirmed Ha, [O 
iii]/H/3, [O n], and Lya lines are summarized in Table 1. 
We note that a small subset (9 NB921 and 10 NB973) of 
our spectroscopic confirmations do not correspond to a 
strong emission line entering our NB filter. These could 
be explain by (1) photometric errors, or (2) an old stel- 
lar population, which would decrease the z'-band flux 
relative to flux at the NB wavelength (see Section 3). 

We illustrate in Figure 3, the A(z'-NB) colors for 
spectroscopically confirmed sources. The NB921 spec- 
troscopic sample spans six magnitudes in luminosity and 
three magnitudes in NB921 excess. Because of sparser 
sampling, the NB973 spectroscopic sample covers four 
magnitudes in luminosity and one magnitude in NB973 
excess. They show that NB921 (NB973) excesses of 0.2 
mag (0.25 mag) are indeed due to the presence of weak- 
but real-emission lines. 

3.1.2. Dual Line Detection with NB Filter Pairs 

In addition to spectroscopy, another independent ap- 
proach showing that the NB selection technique effi- 
ciently identifies line-emitting galaxies, is a second NB 
filter that targets another emission line at the identical 
redshift. Coincidentally, the NB704 filter (centered at 
7046A with a FWHM of 100A) is able to accomplish this 
in conjunction with NB921. Others (e.g., Hippclein et 
al. 2003; Nagao et al. 2008; Nakajima et al. 2012) have 
exploited similar dual-line instances to confirm emission- 
line galaxies at a particular redshift. When Ha (H/3) is 
redshifted into the NB921 filter at z = 0.40 (z = 0.89), 
[O in] ([O ii]) enters the NB704 bandpass. We illus- 
trate in Figure 4 the redshift range probed by these 

http: //www. astro .princeton. edu/~rcool/hsred/. 
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TABLE 1 

Summary of Spectroscopic Samples 



Filter 
(1) 



Targeted 
(2) 



Confirmed 
(3) 



Ambiguous 
(4) 



DEIMOS 

(5) 



MMT 
(6) 



FOCAS 
(7) 



FMOS 2 
(8) 



N(Ba) JV([0 m]) N([0 n]) N(Lya) 



(9) 


(10) 


(11) 


(12) 


81 


80 


104 


21 


19 


12 


32 


1 



NB921 
NB973 



522 
204 



295 
104 



227 
100 



88 (127) 
6 ( 13) 



113 (164) 
61 ( 80) 



36 (48) 
6(9) 



68 (212) 
30 (105) 



Note. — The number of targeted, spcctroscopically confirmed, and unclassified galaxies are given in Cols. (l)-(3). Values in parentheses denote 
the total targeted sample with each instrument. Confirmed galaxies that are neither Ha, [O III], [O II], or Lya are discussed in the text. 

* Note that a subset of the samples were targeted with different instruments. These numbers account for overlapping targets. 

2 The FMOS reduction is not complete, as improvements for optimal spectral stacking and detection of weak emission lines have yet to be implemented. 



Thus, these confirmed numbers represent lower limits. 
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Fig. 4. — Redshift coverage for dual emitters. Dashed (blue) 
and dotted (red) lines show the NB704 and NB921 bandpasscs, 
respectively. The product of the two filters is shown by the solid 
lines. The left (right) figure illustrates the [O m]-H<* ([O irj-H/3) 
combination for z m 0.40 (z ~ 0.89). 

two NB filters for the [O in]-Ha and [O ii]-H/3 com- 
binations. Cross-matching a sample of 1421 3cr-detected 
NB704 excess emitters (L07) with the NB921 emitting 
sample discussed, we found 241 dual NB704-NB921 emit- 
ters. These dual emitters are overlaid on Figure 3 as 
crosses. Among the dual emitters, 80 were spectroscop- 
ically targeted, and only nine of them have unidentified 
redshifts. Among the remaining 71 sources, 66 of them 
(93%) were Ha emitters. There are only two cases (3%) 
of dual [O n] and H/3 at z « 0.89. Since these lines are 
much weaker (by a factor of 5 to 10; see Hicks et al. 2002) 
than [O III] and Ha for typical star-forming galaxies at 
these redshifts, it is not a surprise that the majority of 
spectroscopically confirmed dual emitters are at z ~ 0.4. 

3.2. Identifying [0 II/ Emitters 

Rest-frame optical colors were previously used to dis- 
tinguish the three most common rcdshiftcd emission 
lines, [O II], [O III], and Ha, entering into our narrow- 
band filters (L07). A small sample of spectroscopically 
confirmed sources is used to calibrate and optimize the 
broad-band color selection. With a much larger spectro- 
scopic sample, we adopt an improved hybrid approach to 
select our [O n] samples for our analysis. First, spec-z's 
provide unambiguous classifications of NB excess emit- 
ters. For the NB921 (NB973) sample, 104 (32) candi- 
dates are definite [O n] emission-line galaxies. These 
spectroscopic confirmations represent «10% of the final 
[O n] samples. Next, we apply two-color photometry to 
select [O n] emitters. Illustrated in Figure 5 are the i?c~ 
i' and i'-z' con t colors for our excess emitters, where z' CO nt 
is the "true" continuum at ~9200A after using the NB 
photometric measurements to remove contamination of 



2 ' cont = z' + 2.51og^l + ^4pj, (6) 

where 955A is the FWHM of the z'-band filter. This is 
especially important since [O in] emitters with very large 
EWs will have unusually red i'-z' colors, which can be 
distorted enough to cross our color selection boundaries. 
The median differences in z' and z' cont are 0.09 mag 
(NB921) and 0.14 mag (NB973), with average differences 
of 0.11 mag (NB921) and 0.24 mag (NB973), and can be 
as large as 1.25 mag (NB921) and 1.62 mag (NB973). 
The broad-band color choices were made in L07, and 
were designed to sample both sides of the Balmer/4000A 
break at z ~ 0.9, while for higher redshift, probing blue- 
ward of the break. In this figure, we have already ex- 
cluded Ha emitters with our previous BRqi' two-color 
selection (Ly ct al. 2012). To select [O n] NB921 excess 
emitters with optical colors, we adopt: 

R c -«'<0.05, or 
R c -i'<l.ll(i'-z' cont )-0M. 

These color criteria are similar to those adopted in L07. 
With a larger spectroscopic sample, we determined that 
they only miss 1 of 100 (1%) spectroscopically identified 
[O n] NB921 emitters. Also, only 2 of 70 (3%) identi- 
fied [O in] emitters meet the [O n] color selection. For 
selecting [O n] NB973 excess emitters, we adopt: 

Rc-i'< 0.20, or 

Rc - i' <0 .45(i' -z' COILt )+ 0.22. 

Based on the above broad-band color selection, combined 
with our spectroscopic redshifts, we have identified 933 
[O n] NB921 at z sa 1.47 and 328 NB973 at z k, 1.62. We 
will now test these redshift determinations by comparing 
them with photometric redshifts. 

3.3. Photometric Redshifts for NB Excess Emitters 

While spectroscopy is not available for all galaxies, 
photo-z's can be determined from their individual SEDs. 
In a previous SDF study (Ly et al. 2011), we constructed 
a very large sample of ^200,000 objects in SDF with pho- 
tometric measurements in twenty bands: FUV, NUV, 
UBVR c i'z'z b z r JHK, five narrow-band filters (NB704, 
NB711, NB816, NB921, and NB973), and two intermedi- 
ate bands at ^6500A. By cross-matching our NB excess 
emitter catalogs against the ultra-deep SDF photomet- 
ric catalog, we have full photometric measurements for 
2132 of 2361 NB921 excess emitters and 1040 of 1243 
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Fig. 5. — Rc—i' and i'— z'cont colors for NB921 (left) and NB973 (right) emitters shown as black squares. Ha emitters, identified in Ly 
ct al. (2012), have been excluded from this plot. Spectroscopically confirmed galaxies are shown as triangles with green and blue colors 
representing [O III] and [O Ii] emission entering the NB filter, respectively. 
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Fig. 6. — Photometric redshifts for NB921 (top) and NB973 emitters (bottom) derived from twenty-band photometry. Overlaid are the 
NB921 and NB973 filter profiles at the redshifts where Ho (red), [O III] (green), H/3 (purple), and [O ll] (blue) enter the NB filters. Samples 
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i'—z' color (see Section 3). Vertical dashed lines show the FWHM of the filters. 



NB973 excess emitters. The photo-z's for these sources 
were derived using the Easy and Accurate .Zphot from 
Yale (EAZY; Brammer et al. 2008). We illustrate these 
photo-z's in Figure 6, and compare them against our 
spectroscopic redshifts in Figure 7. There are good agree- 
ments between the photometric and spectroscopic red- 
shifts: the rms of {z phot -z spcc ) / (I + z spcc ) = 1+ ^ ^ is 
0.0039 for NB921 and 0.0083 for NB973. This is not sur- 
prising, because the EAZY fitting routine includes emis- 
sion lines, which are needed to match the observed broad- 
and narrow-band photometry. It appears that the photo- 
metric redshifts derived for Ha and [O n] NB973 excess 
emitters are systematically off by The cause is 

unknown. 



Not including NB filter measurements, the photo-z dis- 
tribution shows three peaks near the expected location 
for Ha, [O in], and [O ii]. However, the photo-z uncer- 
tainties are much larger, for example, [O ii] emitters have 
Zphot's between 1.4 and 1.8. 

We define a Zphot outlier as having 



> 0.15. 



There are 56/250 (22.4%) and 15/91 (16.5%) failures for 
the NB921 and NB973 samples, respectively. As illus- 
trated in Figure 8, the majority of these z p hot failures are 
associated with galaxies with high observed EWs. For ex- 
ample, in the NB921 sample, there are almost two dozen 
spectroscopically confirmed [O ill] emitters at z « 0.84 
with EWobs ^ 300A. All of them have incorrect photo- 
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z's, mostly z p hot ~ 0.4. We determined that the spec- 
tral templates used by EAZY do not allow for extremely 
strong [0 in] emission compared to Ha (i.e., [O in] /Ha 
flux ratio greater than unity). As a result, EAZY finds 
a better fit by assuming that the large NB921 excess is 
due to Ha. As L07 have shown, the [O in] /Ha ratio from 
z 0.4 dual emitters were found to be on average, unity. 
EAZY will need to be revised to account for the recently 
observed fact that star-forming galaxies at high-z have 
very large [O in] /Ha ratios, perhaps due to low gas-phase 
mctallicity (Atek et al. 2011). A similar problem occurs 
for the spectroscopically confirmed Ha emitters where 
they are placed at z ~ 0.1, a redshift which puts Ha in 
the NB704 filter. However, the NB704 excess is in fact 
due to [O in]. This problem occurs when the observed 
Ha+[N n] EW is above lOOA. 

Other studies have examined photo-z's for emission- 
line galaxies and found less catastrophic failures. For ex- 
ample, Xia et al. (2011) found very few catastrophic fail- 
ures between photometric and grism redshifts. However, 
it is not clear what fraction of their sample targeted high 
EWs, so a direct comparison is not possible. We note 
that there are fortunately very few [O n] NB921 failures. 
Only 3 of 36 (8%) have incorrect z p hot- The higher suc- 
cess for the [O il] selection is because, blueward of 3600A, 
the SED is simple: no additional emission lines enter our 
optical bands, and the Balmer/4000A break is shifted 
past 9000A. 

We note that among the photo- z sample, 77% and 58% 
have photo-z's consistent with having [O n] entering the 
NB921 and NB973 filters, respectively. These numbers 
should be taken as lower limits to the real success of the 
[O n] color selection, because many of the apparently 
discrepant photometric redshifts were highly uncertain, 
based on rather limited multi-band data. If we restrict 
our consideration to only the most robustly determined 
photo-z's, we would find an even higher success rate of 



l+z B 



as a function 



100 

EW (Angstrom) 

Fig. 8. — Photometric redshift accuracy, 

of NB921 observed EW. Ha, [O III], and [O o] emitters are denoted 
as red asterisks, green squares, and blue triangles, respectively. A 
significant number of the z = 0.84 [O in] emitters have incorrect 
photo-z's due to high [O Iirj/Ha flux ratios, which are not consid- 
ered in EAZY spectral templates. These galaxies, with high [O III] 
EWs, are believed to have low gas metallicities. 

identifying [O ii]-emitting galaxies. 

4. RESULTS 

4.1. [On] Equivalent Width Distributions 

Figure 9 shows the distributions of the rest-frame [O n] 
EW (EW ), for galaxies detected at > 3<7 in the NB921 
and the NB973 filters. The distribution for NB921 
emitters further emphasizes the excellent sensitivity our 
NB921 imaging provides even for galaxies with relatively 
weak [O il] lines, with EWo ~ 15A. Although our abil- 
ity to detect [O n] falls off for EW 's below —25 A, the 
broad peak of the observed detections around an EWo 
of 43 A suggests that our NB921 survey misses relatively 
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Fig. 9.— Rest-frame [O II] EW distributions for our NB921 (solid 
line) and NB973 (dashed line) samples. The dotted lines indicate 
the EW selection limits of our surveys. Also shown are the distri- 
butions when the samples are limited to NB excesses of > 4cr (red) 
and > 5<r (blue), analogous to Figure 6. 

few galaxies, since most of them have EWs larger than 
our sensitivity limit. One piece of evidence that supports 
this is the comparison of the EW distributions with more 
conservative NB excesses limits (4cr and 5a; see Figure 9). 
In these cases, the distributions still have the same shape: 
peaking at the same EW with the lower EW slope un- 
changed. If there was a significant number of 15-25A 
[O n] emitters detectable at 3-5cr, we would see a more 
rapid decline prior to our EW selection limits. 

Recently, Kornei et al. (2012) also investigated the [O 
n] EW of star-forming galaxies at somewhat lower red- 
shifts of z — 0.7-1.3. They find a similar median [O n] 
EWo of 45A, and also find a similar decline towards lower 
EWs. This result provides further independent support 
for our claim that our NB921 survey misses few galaxies 
in the redshift range for [O 11] , because that line is nearly 
always present. 

In addition, another NB survey that examined the Ha 
EWs at z = 2.2 finds that the Ha EWo distribution peaks 
at 140A, «120A above the survey limit, and intrinsically 
falls off rapidly towards lower EWs (Lee et al. 2012). 

4.2. Comparison of [O ll] and UV Continuum 
Luminosities 

For galaxies at any appreciable redshifts (z > 0.5), the 
most common observable used to estimate their current 
star formation rates (SFRs) is their UV continuum emis- 
sion. This works because in all but the oldest stellar pop- 
ulations, the far-UV continuum (1500A rest wavelength) 
should be dominated by the hottest, and thus shortest- 
lived, main-sequence stars (O and B types). Since the 
hottest of these produce most of the ionizing photons in 
a galaxy, it is expected that as the rate of current star 
formation increases, so docs the UV continuum as well as 
the emission line luminosity from H II regions. The "gold 
standard" measure used to count ionizing photons is the 
Ha recombination line, but this is redshifted beyond the 



sensitivity of the best CCD spectrographs at z ~ 0.5. 
To study higher redshifts, secondary emission lines from 
ionized gas have been explored at shorter wavelengths. 
The most "popular back-up" indicator of SFR is the lu- 
minosity of the [O 11] doublet. 

Figure 10 shows the observed luminosities of the [O 
ll] emission line as a function of the observed far-UV 
continuum at 1500 A. The solid line shows the least- 
squares regression after binning our data in £(1500): 
logL(On) = 16.23 + (0.89 ± 0.22) log £(1500). Here, 
£(1500) is determined directly from the [/-band photom- 
etry 



log 



£(1500) 



erg/s/Hz 



(Mu + 48.6) 
^5 



40.08, 



(7) 



where My = [7+2.5 log (1 + z)-5 log (d L /10 pc), and d L 
is the luminosity distance, 1.06 x 10 4 Mpc for z = 1.47. 
And the [O n] luminosity is determined from narrow- 
band and broad-band measurements: 



£(On) = ANB 



A,NB 



fx 



1 - (ANB/Az) 



(8) 



In fitting the £(1500)-£([O ll]) relation, we have ex- 
cluded measurements below log £(1500) = 27.75, since 
the [O n] measurements are biased against galaxies that 
are fainter than the current [O n] depth permitted by 
our NB imaging. 

Both of these observables have been calibrated to pre- 
dict the SFR of a galaxy (Kennicutt 1998), if the quan- 
tities can be corrected for interstellar extinction. Kenni- 
cutt (1998)'s SFR relation for £(On) is based on the cal- 
ibration for Ha. In his sample of nearby spiral galaxies, 
he estimated £(On) = 0.45 £(Ha). 12 He recommends an 
average Ha extinction correction factor of 2.5 (Kennicutt 
1992). Thus, we have 



SFR = 2.0 x 10" 41 £(On) 



L(Ha) 



£(On) 



(9) 



For a Salpeter (1955) (hereafter Salpetcr) initial mass 
function (IMF) , the far-UV continuum and SFR are also 
related by a simple linear proportionality: SFR = 1.4 x 
10~ 28 £(1500) x io°- 4A ( 150 °), where ,4(1500) = 1.37 mag- 
assuming a Calzetti et al. (2000) dust curve with j4(Ha) 
= 1 mag and E(B - V) star = 0.UE(B - V) gas; If we 
require that both observable diagnostics must yield the 
same SFRs, then we would expect a simple linear relation 
between observed £(Oil) and £(1500): 



£(On) = 2.5 x 10 id £(1500) 



£(Oii) 
L(Ha) 



(10) 



This linear relation is shown in Figure 10 for two values 
of the observed [O n]/Ha flux ratio, 0.5 and 1.0. This 
comparison suggests that, if the Kennicutt (1998)'s SFR 
relations are correct for £(1500) and for £(Ha), then 
our data indicate that the [O n]/Ha flux ratio is typi- 
cally 0.75 for the majority of our sample. Sobral et al. 
(2012) also compared the [O ll] and Ha flux ratio at the 
same redshift and found a typical ratio of 0.4-0.8. For 
a subset of our sample with high UV luminosities, the 

12 L(On) and L(Ho) are observed (i.e., not de-reddened). 
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Fig. 10. — Observed [O Ii] and UV luminosities of z = 1.47 [O n] emitters. The greyscale represents the density of points in each bin. 
Purple filled circles shown the median at each UV luminosity bin with error bars representing 68% of the [O Ii] measurements. Overlaid 
as dashed lines are the cases where the observed [O Il]/Ha flux ratio is 0.5 and 1.0 following Equation 10. The histograms collapse the 
distributions into one dimension with blue and red colors indicating those [O Ii] emitters meeting the BX/BM and sBzK selections (discussed 
further in Sections 4.5.1—4.5.2), respectively. 



[O n]/Ha ratio decreases towards 0.5. This is consistent 
with the same decline in [O ii]/Ha flux ratio with blue 
stellar luminosity found by Jansen et al. (2001) in local 
galaxies. Illustrated in Figure 10, the scatter of the [O 
ii] luminosity relative to the UV luminosity suggests that 
[O Ii] can be used as a SFR indicator with a reliability 
of a = 0.23 dex. 

4.3. Average Spectral Energy Distributions of [O 
il]-Emitting Galaxies 

We have combined our multi-wavelength photometry 
from hundreds of [O n]-emitting galaxies to construct 
representative broad-band SEDs. We limit our results 
to the NB921 sample, since the results are better con- 
strained with the larger and deeper sample. Our filters, 
U through IRAC 1 and 2 (Fazio et al. 2004), span rest- 
frame 1500A-1.8/^m, more than ten times the spectral 



coverage that was possible with the Suprime-Cam data 
alone, that was used by L07. 

In Figure 11 we show the median SEDs of z = 1.47 
[O ii] emitters, dividing the sample between high- and 
low-equivalent width of [O ii] emission (greater or less 
than 106A) and between low and high [O ii] flux (above 
and below 1.49 x 10~ 17 erg s _1 cm -2 ). For comparative 
purposes, all SEDs are normalized to 23.0 mag in the i' 
band. The error bars show the range encompassing 68% 
of the data in each waveband. 

The average U-to-i' photometry (corresponding 
roughly to rest wavelengths between 1500 and 3300 A) 
can be approximated by a power law (F u cx v a ) with 
a = —0.64, except for a marginal but systematic dip in 
the intermediate 6000 A filter, which might possibly be 
attributed to the broad 2200 A peak in the interstellar 
reddening law. On the red side of the i' band, the SEDs 
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rise much more rapidly, due to the Balmer and Ca II H 
and K breaks with an average i' — K color of 1.2 mag. 
Redward of the jumps, the rest-frame optical continuum 
is relatively flat. This is supported by the IRAC [3.6] 
and [4.5] photometry. We note also that at this redshift 
of z = 1.47, [O III] and H/3 enter the J-band, and thus 
can alter these colors if these emission-line EWs are ex- 
tremely high, as seen in other surveys (e.g., Atek et al. 
2011). Indeed, if we compare the average J flux of the 
high-EW galaxies to a typical model SED going through 
the surrounding filters, it appears to be a few tenths of a 
magnitude brighter. This is attributable to the contam- 
ination from the [O in] and H/3 lines. 

As found by L07, the average SED of the high-EW sub- 
sample is systematically bluer than of the moderatc-EW 
subset: the observed < U — i' > and < i' — K > colors 
are 0.43±0.27 and 1.13±0.35 and 0.58±0.32 and 1.28± 
0.50 for the high- and low-EW sub-samples, respectively. 
This is reasonable, since the former galaxies, with specific 
SFRs several times higher, should have a larger fraction 
of younger, blue stars. We also find that galaxies with 
higher [O n] emission-line fluxes are found to be redder 
than those with weaker [O n] flux. We argue that this is 
due to a selection effect that more massive galaxies tend 
to be redder will have higher SFRs and thus higher [O II] 
fluxes. We next quantify these trends with fits of model 
stellar populations to the SEDs. 

4.4. Stellar Populations from SED Modeling 

We modeled the fifteen-band SEDs of our [O n] emit- 
ters to derive stellar masses, SFRs, dust reddening 
(Ay's), and stellar ages. At z ~ 1.5, the data span 
rest-frame 600A-8800A. The spectral synthesis models 
are discussed further in Ly et al. (2011). In brief, we use 
the IDL-based code, Fitting and Assessment of Synthetic 
Templates (FAST; Kriek ct al. 2009), with exponentially 
declining r star-formation history models from Bruzual 
& Chariot (2003). We adopt a Salpcter IMF, and dust 
attenuation follows the Calzetti et al. (2000) reddening 
formalism. A Chabrier (2003) IMF was also considered; 
however, the resulting x 2 fits (ages, r's, and Ay's) were 
no different from those determined using a Salpeter IMF. 
The only difference was a lower stellar mass and SFR 
by 0.25 dex (a factor of 1.8), which is expected due to 
the smaller proportion of low-mass stars. Therefore, we 
adopt the Salpeter IMF for the remainder of this paper. 
To correct for neutral hydrogen intergalactic medium ab- 
sorption, we follow the prescriptions described in Madau 
(1995). 

The grid of models consists of log(r/yr) = 7.0, 8.0, 
9.0, and 10.0. These r values were selected to in- 
clude a bursty, intermediate, and roughly constant star- 
formation history. In addition to the star-formation his- 
tory, the grid of models spans a range of log(age/yr), 
between 6.0 and up to the age of the universe at a given 
redshift (in increments of 0.1 dex), and dust extinction 
with Ay = 0.0-3.0 mag with 0.1 mag increments. We 
only consider solar metallicity models, since distinguish- 
ing different metallicities is nearly impossible with broad- 
band SED fitting (see e.g., Pforr et al. 2012). 

In these fits, the rcdshifts are accurately fixed, cither 
from the spectroscopic values when available, or from 
the presence of [O 11] in the NB bandpass (z = 1.467 



or z = 1.617). Since these spectral synthesis models do 
not include nebular emission lines, we correct the z'-band 
measured fluxes to remove the [O 11] emission contribu- 
tion to these filters. We note that by having the pho- 
tometry as is, the stellar masses were higher by 0.2 dex, 
the stellar ages were older by 0.1 dex, and the SFRs were 
lower by 0.3 dex. 

In Figure 12, we illustrate the stellar masses and ages 
for 863 z « 1.47 [O 11] emitters. We find that these [O 
11] emitters span 2 dex in stellar age with a median of 
^100 Myr and 3 dex in stellar mass with a median of 
3 x 10 9 M . In modeling the broad-band SED, there is a 
known degeneracy between the fitted stellar age and dust 
reddening in describing the overall spectral slope, partic- 
ularly in the UV. In particular, one can make galaxies 
appear younger but with significant dust and vice versa. 
We caution that fits with very young stellar populations 
(<10 Myr) may not be reliable because they also give 
very large extinction (Ay = 1-3 mag). 

We find that the galaxies with observed [O 11] EW 
above 106A have on average a factor of two less stel- 
lar mass than the low-EW galaxies. In addition, galaxies 
with [O 11] fluxes above 1.49 x 10~ 17 erg s _1 cm~ 2 are 
more massive by a factor of 5. These results are consis- 
tent with the average SED results discussed in Section 4.3 
since more (less) massive galaxies are typically expected 
to have redder (bluer) colors. 



4.5. Testing Color- Selection Techniques at z > 1 

The two techniques commonly used to select star- 
forming galaxies at redshifts of 1-3 are the "BX/BM" 
(Adelberger et al. 2004) and "sBzK" (Franx et al. 2003; 
Daddi et al. 2004) methods. The former approach, anal- 
ogous to the Lyman break method (Steidel & Hamilton 
1992), finds young galaxies with strong star formation 
and moderately red U n — G color due to the Lyman con- 
tinuum break beginning to enter into the U n filter. The 
latter method is intended to identify z > 1.4 galaxies 
when the Balmer/4000A break occurs between z' and K. 
Because of the need for A"-band data, the BzK method 
is generally biased towards more massive galaxies. These 
methods were developed to help "fill" the historically 
known "redshift desert" where wide-field deep UV imag- 
ing to identify z ~ 1.5 Lyman break galaxies was not 
available until more recently (Ly et al. 2009; Hathi et al. 
2010; Oesch et al. 2010; Basu-Zych et al. 2011; Haberzettl 
et al. 2012). These color selections and their respective 
samples identified for the SDF are discussed in Ly ct 
al. (2011), where they illustrated that a wide range of 
the stellar populations of galaxies at z > 1 is probed by 
combining these two complementary methods. 

With deep multi-wavelength coverage for the SDF, we 
examine which of our [O 11] emitters are classified as UV- 
selected "BX/BM" galaxies and/or NIR-selected BzK 
galaxies. 

Since we argued above that most galaxies at z = 1.47 
have detectable [O 11] emission lines, we can use our large 
NB samples to test the completeness of the two main 
broad-band color search methods available at the same 
redshift. We find that both suffer from some incomplete- 
ness, which is mostly eliminated by combining the two 
methods. 
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(right). Each SED is constructed from roughly 350 galaxies where the sub-samples are constructed above and below the median EW and 
[O n] flux. Error bars represent 68% of the sub-samples, and blue and red colors denote the samples with bluer and redder rest-frame UV 
colors. 




niverse ot z= 1 .47 



□ □□□ 



□ □□ 



] □□ 
]□ □ 

: □ 



□ 

□□ 

□ □ □□□[_ 
□□ □ □□ 
□ □□□ 
□ □ □□□ 

□□□□□□□□□□□□□L 
□□ 

□ □□□□□□[ ~ 



□ □ 
□□ 

_D □ _ 
□ □ C 
WW 
□ 



□□ □ □□□□□ □ 
□□□□□□□□□ □ r 
□□□□□□□□□□□□□□ 



□ □ 
□□ □ 
□ 



□□□ 
□ 



□ □ 
□□ □ 
□ □ □ 



□ □ 

□ 

□□ □ □ 




50 100 150 



log(M[M,„,]) 




□ □ 




□ □ □ 



□ m 

■ m 

□ □ □ □ cm 



an cd a 
m m m 
c m □ □ □ □ cd c 
CD mm rm □ rm 
m □ m □ □ 




logtMfM^]) 



Fig. 12. — Stellar ages, stellar masses, and [O II] observed EWs for the z = 1.47 [O II] NB921 emitters. These quantities were derived 
from SED modeled fits that assumed a Salpeter IMF, an exponentially declining star-formation history, and Calzctti et al. (2000) internal 
dust reddening. The histograms collapse the distributions into one dimension with blue and red colors indicating those [O II] emitters 
meeting the BX/BM and sBzK selections, respectively. 



4.5.1. The BX/BM Sample 

Figure 13(a) shows the U n — G and G — 1Z colors for 
both of [O n] emitting samples. Since the SDF does not 
have these exact filters, we derive these colors by using a 
combination of the U BVRci' filters that are measured. 
The color transformation was determined by using stellar 
population synthesis models representative of UV lumi- 
nous galaxies (Ly et al. 2011), and consists of: 



U n = U, 



(11) 



G = BV 



-2.5 log 



K = R c i' = -2.5 log 



xi/b + (1 - xi)fy 

3630 ,uJy 
x-iSr + (1 - x 2 )Jj 



3630 /LtJy 



and (12) 



(13) 



Here, fx is the flux density per unit frequency (erg s 
cm~ 2 Hz" 1 ) in band "X", and X\ = 0.314 and x 2 = 
0.207, as discussed in Ly et al. (2011). 

Limiting our consideration to i?c ^ 25.5, the typical 
depth adopted by past and current BX/BM surveys, we 
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Fig. 13. — (a) U - BV and BV - Rci' colors for [O II] emitters. Red symbols denote sources fainter than Rq = 26.0. BM, BX, and 
[/-dropout galaxies are selected by the three regions defined by the solid lines (from bottom to top, respectively). Most galaxies are 
classified as BM, which is expected for z = 1.47 and 1.62. (b) B—z' and z'-K colors for [O II] emitters. Red symbols identify sources that 
are undetected in the ii"-band at the 3cr level (upper limits on z'-K). Star- forming BzK galaxies are classified as those lying above the 
diagonal line. Both panels show blue crosses to denote the galaxy population selected by the other method. The NB921 and NB973 [O II] 
samples arc distinguished by squares and triangles, respectively. Spectroscopically confirmed [O II] emitters are shown as orange circles. 



find for the full z = 1.47 [O n] sample that 47% are BX 
or BM. Likewise, 82% of the z = 1.62 [O n] emitters are 
either BX or BM. We summarize these statistics, as well 
as those for R c < 26.0, in Table 2. 

4.5.2. The sBzK Sample 

The B — z 1 and z' — K colors for the [O n] emitters are 
illustrated in Figure 13(b). Since the SDF is not fully 
covered (80%) with UKIRT//v-band observations to a 
3cr sensitive depth of 23.9 mag, we limit the sample of 
[O n] emitters to those that arc in the sensitive region. 
Among the 863 z = 1.47 (313 z = 1.61) [O n] emitters, 
750 (264) sources are located within the sensitive region. 
Requiring a 3cr detection in K, the samples are further 
limited to 402 (NB921) and 174 (NB973). Daddi et al. 
(2004) defines star-forming BzK (sBzK) as galaxies with 



K >(B 



0.2. This yields 84% of [O n] NB921 



and 91% of [O n] NB973 emitters as sBzK galaxies. None 
of these [O n] emitters have colors that would classify 
them as passive BzK. Thus both color methods, espe- 
cially BzK, work very well to find emission-line galaxies. 

We emphasize that 80 of 576 (14%) .ftT-band detected 
[O n] emitters have blue z 1 — K colors preventing them 
from being selected as sBzKs. As illustrated in Fig- 
ure 13(b), a large fraction of these has been spectroscop- 
ically confirmed to be at z = 1.47 or z = 1.62. Lee et 
al. (2012) have also shown that for a sample of z = 2.2 
Ha emitters, 8% are also missed because of their blue 
z' — K colors. They argued that the sBzK selection 
fails to identify galaxies with relatively younger stellar 
population due to little or no Balmer break. Barger et 
al. (2008) also examined the BzK selection with a large 
spectroscopic sample and find that the sBzK selection 
is effective at identifying high redshift galaxies, but may 
suffer from significant (33%) contamination from z < 1.4 



interlopers. 

4.5.3. Stellar Population Characteristics of Color-Selected 

Galaxies 

To illustrate the diversity of the [O n] sample and the 
selection biases associated with the BX/BM and sBzK 
techniques, we plot their distributions of stellar masses, 
ages, and [O n] EW in Figure 12. It is apparent that the 
sBzK selection (shown in red) is well designed to identify 
almost all of the massive [O n] emitters, but becomes 
seriously incomplete below 3 x 10 9 M Q . This incom- 
pleteness is caused by the current depth in the .ftT-band. 
By contrast, the BX/BM selection (shown in blue), suc- 
ceeds at identifying most of the low-mass galaxies, but 
misses about half of the massive population. In addi- 
tion, we present the first determination of [O n] EW 
distributions for large representative samples of BX/BM 
and sBzK galaxies. The BX/BM sample is skewed to- 
wards higher observed EWs (average of 105A) compared 
to sBzKs (average of 91 A). 

We illustrate in Figure 14(a) the observed [O n] 
emission-line EW as a function of Rq-z' . It is clear that 
the BX/BM and sBzK techniques are complementary by 
identifying the bluest and reddest galaxies, respectively. 
Roughly 25% of the [O n] emitters arc identified using 
both color selections, while 34% and 20% are identified 
by only the BX/BM and sBzK selections, respectively. 
Figure 14(a) also shows that -20% (154/750) of the [O n] 
sample is not identified by either method. These galax- 
ies are generally of lower luminosity where 43% of them 
are fainter than Rq =26.0 mag and 84% arc undetected 
(at 3a) in K. These galaxies have typical limits on their 
stellar mass of — 3 x 10 8 M Q . We note that the [O n] 
emitters galaxies with Rc < 26.0 but which still miss the 
BX/BM selection have typical stellar masses of 2 x 10 9 
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TABLE 2 

Summary of BX/BM and sBzK Results for [O ii] Emitters 



Filter 


N 


-Nphot 


N(R < 25.5) 


N(R < 26.0) 


UV(R<25.5) 


UV(R<26.0) 


jV(K-det) 


iV(sBzK) 


JV(Either) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


NB921 


933 


863 92.5% 


585 67.8% 


741 85.9% 


406 47.0% 


520 60.3% 










(814) 


(750 92.1%) 


(505 67.3%) 


(641 85.5%) 


(348 46.4%) 


(446 59.5%) 


(402 53.6%) 


(338 45.1%) 


(596 79.5%) 


NB973 


328 


313 95.4% 


286 91.4% 


307 98.1% 


255 81.5% 


272 86.9% 










(278) 


(264 95.0%) 


(238 90.2%) 


(258 97.7%) 


(210 79.5%) 


(226 85.6%) 


(174 65.9%) 


(158 59.8%) 


(248 93.9%) 



Note. — Values reported in parentheses refer those sources that fall in the sensitive iC-band region ("K-det"). Except for those reported in 
Col. (3), all percentages are relative to -/V p h ot . In Cols.(6)-(7), "UV" refers to the combination of "BX" and "BM" selections. 



(a) 




N = 750 
only) = 258 
® NfsBzK only) = 150 
N(both) = 188 
O [Oil] spec-z 
® ■ 



,3 .0. 



■'f.<®^<. ' ■' ' 

■ 4 ®; ®f. ■ 1 ■ .■ 




0.0 0.5 1.0 1.5 41 -° 41 - 5 42 -0 42.5 

R c - z' log(L oss ([0 ll])/erg s') 

Fig. 14. — (a) [O Ii] observed EW as a function of Rc~z' color for NB921 emitters. [O n] emitters that are selected purely as "BX/BM" 
{Rq < 26.0 mag) and sBzK are shown as blue and red squares, respectively. Those identified using both color selection techniques are 
shown as green squares, and those that do not meet both selections are shown as black squares. Spectroscopically confirmed [O Ii] emitters 
are shown as black circles, (b) The probability distribution of [O Ii] observed luminosity, as shown by the solid black line. The red dashed 
line illustrates the distribution for [O Ii] sBzK galaxies. Likewise, the blue dotted line shows the [O Ii] BX/BM galaxies. The combination 
of both method is shown by the shaded histogram. Percentages are reported relative to the full [O Ii] sample. 



M©, which is 0.25 dex more (0.86 less) massive than the 
BX/BM (sBzK). This suggests that a subset of inter- 
mediate mass star-forming galaxies are missed by color 
selections. 

Finally the distribution of [O n] luminosity for the 
sBzK and BX/BM samples are shown in Figure 14(b) 
where it illustrates that either of these color selections 
will miss 40%-50% of the [O n] sample. Even combining 
both methods, roughly 20% of the faintest [O n] emitters 
are missed. 

5. CONCLUSIONS 

Using deep NB imaging with Suprime-Cam, we identi- 
fied ^1,300 [O ii] emitting galaxies at z=1.5-1.6. Follow- 
up optical and near-infrared spectroscopy confirmed that 
the two-color selection of [O n] emitters is highly reliable: 
Only 1% of NB921 [O n] emitters appeared to be missed 
while 3% of z = 0.84 [O in] emitters contaminate the [O 
ii] NB921 color selection. Our survey has investigated 
the properties of [O n]-emitting galaxies. We find: 

1. The average rest- frame [O Ii] EW of our z = 1.47 
sample is 45A, consistent with another study at 



z > 0.7. The EW distribution declines gradually 
from the peak towards lower EWs, suggesting that 
very few z ~ 1.5 galaxies lack [O ii] emission. 
We compared the observed [O n] and UV luminosi- 
ties and find a strong correlation with log L(0 ii) oc 
(0.89 ± 0.22) log 1,(1500). The observed scat- 
ter (0.23 dex) of individual [O n] measurements 
against the UV luminosity suggests that the [O n] 
luminosity can be considered as a coarse SFR indi- 
cator. 

We constructed multi-band SEDs (rest-frame 
1500A-1.8/Ltm), and find that higher EW galaxies 
have slightly bluer colors compared to lower EW 
galaxies. 

We fitted the SEDs with exponentially declining 
"t" models, and find a median stellar ages of 100 
Myr with 13% of the sample older than 500 Myr. 
The stellar masses are typically 10 9 5±0 - 6 M Q . 
We find that [O Ii] emission is present in the full 
census of galaxies, from massive red systems (found 
by sBzK) to smaller, star-bursting systems (found 
by BX/BM). With the sensitivity that we achieved 
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in the NB921 filter, we find that 154 («20%) 
z = 1.47 [O n] emitters are not identified with BzK 
or BX/BM because they are too faint (i.e., stellar 
mass limits of ~ 3 x 10 8 M Q ). A subset of these un- 
identified galaxies that arc bright enough to satisfy 
the BX/BM magnitude limit (R c = 26.0) have in- 
termediate stellar masses (2 x 10 9 Mq) compared 
to the BX/BM and sBzK samples, and suggest 
that some intermediate mass galaxies are currently 
missed by commonly used high-z color selections. 



tional help using the MMT/Hcctospec hsred reduction 
pipeline, and Ryosuke Goto for assistance in the FMOS 
observations. The authors wish to recognize and ac- 
knowledge the very significant cultural role that the sum- 
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